Theoretical and experimental analysis of H2 binding in a prototype metal organic 

framework material 



00 
O 

o 

(N 

Q 
in 



C/3 



o 



> 
o 
o 
l> 

(N 
00 

o 



X 



Lingzhu Kong,-'^ Valentino R. Cooper,^' ^ Nour Nijem,'^ Kunhao Li,"* Jing Li,"' Yves J. Chabal,'^ and David C. Langreth-' 

'Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854-8019 
^Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-611^ 

'^Department of Materials Science and Engineering, 
University of Texas at Dallas, Richardson, Texas 75080 

'^Department of Chemistry and Chemical Biology, 
Rutgers University, Piscataway, New Jersey 08854-8087 

Hydrogen adsorption by the metal organic framework (MOF) structure Zn2(BDC)2(TED) is in- 
vestigated using a combination of experimental and theoretical methods. By use of the nonempirical 
van der Waals density-functional (vdW-DF) approach, it is found that the locus of deepest H2 bind- 
ing positions lies within two types of narrow channel. The energies of the most stable binding sites, 
as well as the number of such binding sites, are consistent with the values obtained from experimen- 
tal adsorption isotherms and heat of adsorption data. Calculations of the shift of the H-H stretch 
frequency when adsorbed in the MOF give a value of approximately —30 cm~^ at the strongest 
binding point in each of the two channels. Ambient temperature infrared absorption spectroscopy 
measurements give a hydrogen peak centered at 4120 cm~^, implying a shift consistent with the 
theoretical calculations. 



Hydrogen storage technology is one bottleneck for the 
utilization of hydrogen as an energy source for mobile ap- 
plications. Metal-organic frameworks (MOFs) comprise 
a rather new class of porous materials in which metal ions 
or clusters are linked by organic units [if. The scaffold 
structure and large apparent surface area make these ma- 
terials potential candidates for hydrogen storage applica- 
tions. One such material, MOF-5, was shown to adsorb 
1.3 wt.-% of hydrogen at 77 K and 1 atm A recent 
effort with Zn2(BDC)2(TED) reported an uptake of 2.1 
wt.-% 0]. Here BDC is benzenedicarboxylate (C8H4O4) 
and TED is triethylenediamine (C6H12N2). To rationally 
design or modify MOF structures to meet storage needs, 
it is critical to understand the interaction between the 
MOF matrix and the adsorbed hydrogen molecules. Fur- 
thermore, the quest to understand the physisorption and 
possible dissociative chemisorption of H2 within a MOF is 
a fascinating problem in its own right. In the present let- 
ter, we show how a relatively recently developed nonem- 
pirical theoretical method 0, Q and infrared absorption 
measurements [gJ can be combined with isotherm and 
heat of adsorption data [3|] to obtain a detailed picture of 
the adsorption of H2 in MOFs. For this purpose we use 
Zn2(BDC)2(TED) f| as a prototype. 

Spectroscopic and theoretical studies have confirmed 
that H2 is weakly bound in MOFs and the binding is 
mainly from long-range dispersive interactions The 
nonlocal electron correlation in this type of interaction 
makes accurate and efficient electronic structure calcu- 
lations very difficult. Quantum chemical methods that 
account for these interactions typically scale poorly with 
system size and are limited to only fragments of the true 
MOF structure. On the other hand, ordinary density 
functional theory (DFT) methods, which are efficient and 
scale well with system size, fail to reproduce the correct 



behavior of the van der Waals interactions, which are im- 
portant in these systems [1]. For example, an H2 bind- 
ing energy of 21.7 meV was obtained for the main H2 
adsorption sites in MOF-5 by use of the generalized gra- 
dient approximation (GGA) within DFT while the 
experimental value is 49-54 meV [l^ ■ 

An alternative, a van der Waals density functional 
(vdW-DF), has been developed by our group and col- 
laborators [1, . It incorporates the van der Waals in- 
teraction into a fully nonlocal and nonempirical density 
functional for the correlation energy, which retains ordi- 
nary DFT's good description of covalent bonding. The 
method has been successfully applied to various sparse 
systems 0, 0, [lH . Here we begin by benchmarking vdW- 
DF on the binding of H2 to an isolated BDC linker, 
where there exist calculations via well established quan- 
tum chemical methods 12, [l3|- It is found that H2 
binds most strongly when pointed along the symmetri- 
cally placed axis through the center of the benzene ring, 
and sitting a distance a little over 3 A away with a bind- 
ing energy of approximately 4 kJ/mol. We find that while 
a standard GGA calculation using the PBE functional 
14| fails to predict more than a third of the interaction 
energy, the vdW-DF [3l not only gives the correct bind- 
ing energy, but quite well reproduces the whole bind- 
ing curve [l3| and the potential energy landscape [l3| , 
though, as is typical [1, 0, |Tl[, the intermonomer equi- 
librium distance is slightly too large. In the full MOF, 
we find that this site's binding energy increases by 3 
kJ /mol due to long range interactions with other parts of 
the MOF — a very nonlocal effect that cannot be repro- 
duced by standard density functionals. More spectacu- 
larly this binding point of the fragment becomes only a 
saddle point in the potential energy surface of the MOF, 
showing that calculations on fragments alone do not nec- 



essarily give relevant information about the binding sites 
in the full MOF, or even useful information as to where 
to look for them. 




FIG. 1: A perspective top view of the unit cell of the MOF. 
The BDC linkers and the metal atoms are evident in this 
view, while the TED pillars are more clearly seen in the two 
subsequent figures. Atomic color scheme: black C, green N, 
red O, blue-green Zn; H atoms are not shown. The locations 
of Channels A (thin blue) and B (thick brown) are schemati- 
cally indicated. The dashed circles indicate places where the 
respective channels continue into the adjoining unit cell. 

The basic MOF structure was taken from experiment 
Based on the single crystal data of the guest-free 
structure the C in TED and O in BDC have a 4- 
fold and 2-fold disorder, respectively. A simplified struc- 
ture eliminating the disorder is used in the actual calcu- 
lations [iBl, as shown in Fig. [TJ The underlying Bra- 
vais lattice is tetragonal with a = b = 10.93 A and 
c = 9.61 A, with a basis consisting of one formula unit of 
Zn2(BDC)2(TED). Our vdW-DF calculations found that 
the strong H2 binding regions are located in two types of 
channel (see Fig. [T]). One such channel type runs in the c 
direction through the length of the crystal, with four such 
channels entering and leaving each unit cell. We denote 
these as Channels A. An energy contour map of one of 
them is shown in Fig. [51 These four channels are not quite 
equivalent in the calculated model because of the differ- 
ing faces presented by the nearest TED and oxygen twist, 
but the differences are relatively small (~1 kJ/mol). A 
second type of channel, which we refer to as Channel B, 
runs in the a or & direction. These run only from one unit 
cell (as we draw it here) to the next, because they would 
otherwise run through the center of the cell, which is a 
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FIG. 2: A diagonal side view of half the unit cell, showing a 
contour energy map of one of the four instances of Channel 
A. The H2 molecule was aligned in the z direction. The units 
of the color scale are kj/mol. 

weak binding region. One of the halves of Channel B is 
illustrated in Fig. [3l There are four halves of Channel B 
per unit cell, i.e. two Channel B's per cell. What is clear 
from the the maps in the Figs. [2] and [31 is that the mean 
binding energy in the maximum regions is 10 kJ/mol in 
round numbers for each of the channels, with Channel 
B perhaps a little higher. Although the maximum bind- 
ing energies in Channels B are affected only very slightly 
(^0.2 kJ/mol) by the orientation of the nearest TED and 
oxygen components, the widths of these channels in the 
xy plane (not illustrated) are more strongly affected. It 
is interesting that unlike what is ap parently true for some 
isoreticular MOFs (IRMOFs) [13], there is little binding 
specifically associated with the metal atoms or in their 
vicinity. 

There are two caveats with respect to the predicted 
binding of ~10 kJ/mol. First, the binding is dependent 
on the H2 orientation (see Table [l|, with a variation of 
~1 kJ/mol. Second, the channels are narrow, and H2 is a 
quantum mechanical object. Its ground state energy will 
be higher than the potential minimum. A full treatment 
is beyond the scope of this letter, but harmonic zero point 
energy estimates indicate that this effect could reduce the 
effective binding by 1 to 2 kJ/mol per II2. 

The above results are consistent with the experimen- 
tal II2 uptake curves at T = 77 K and 87 K [3|], which 
we found to be Langmuir isotherms to a high degree of 
accuracy. The Langmuir isotherm is simply a plot of a 
Fermi function against exp(/x/fcr) {fi is the chemical po- 
tential), combined with the realization that exp(fi/kT) 
is proportional to the pressure P outside the MOF. The 
fact that a single Langmuir isotherm (rather than a linear 
combination of them) fits the uptake data implies that 
the important adsorption sites all have about the same 
binding energy, just as we find. The experimental bind- 
ing energy can be determined both from the isotherm fit 
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and independent heat of adsorption measurements 
each give a value of ~7 kJ/mol. When the quantum cor- 
rections mentioned above as well as thermal excitations 
involving them are accurately calculated and applied, we 
expect that our well-bottom value of ~10 kJ/mol will be 
significantly reduced and closer to experiment. Finally, 
our fit to the Langmuir isotherms gives ^13 adsorption 
sites per unit cell. An inspection of Figs. [2] and [3] implies 
a total of 12, provided that we assign two per unit cell to 
each of the 4 Channels A per cell, and 1 to each of the 
Channel B halves. More importantly, the magnitude of 
the interaction energy increases approximately linearly 
as we load hydrogen molecules on these sites one by one. 

We now turn to the calculation of the H2 stretch fre- 
quencies. We carried out a series of calculations varying 
the bond length of H2, with the center of H2 and the 
host atoms fixed at their equilibrium positions. The re- 
sulting potential energy curve for each configuration was 
used in the Schrodinger equation to obtain the eigenval- 
ues and excitation frequencies [3| • At each of the two 
minimum positions, we performed the calculations for 
all three molecular orientations along x, j/, and z. For 
the Channel A minimum, we calculated the stretching 
frequency along the body-center diagonal as well. We 
express our results in terms of the shifts from our cal- 
culated value of the frequency for free H2 of 4160 cm^^. 
Although the latter is accidentally more accurate than 
expected for the level of approximation used, no use is 
made of this fact, and we consider only the frequency 
shift upon adsorption to be relevant for the comparison 
with experiment. 

The vibrational frequency of the hydrogen molecule 
shifts downward when absorbed into the MOF matrix, 




FIG. 3: A side view of half the unit cell, showing one of the 
four instances of Channel B, each of which extend an equal 
amount into the respective neighboring unit cell. The contour 
map color scale and the H2 orientation are the same as for 

Fig.ia 



as shown in Table ID The relative change from the free 
H2 value, is mainly due to the van der Waals interaction, 
which is weak. The shifts predicted by an application of 
the GGA were found to be substantially smaller still. The 
shifts that we predicted with the vdW functional were 
—28 cm~^ and —25 cm^^ for the preferable binding orien- 
tations in Channels A and B respectively and are smaller 
in magnitude for other molecular orientations. The in- 
teraction energies for different orientations are within 
~1 kJ/mol for each channel, which indicates that all the 
orientations will be well populated at room temperature. 
Furthermore, the translational states, which have a range 
of frequencies of up to ~100 cm"-'^, will be well occupied, 
implying gain and loss contributions of comparable am- 
plitudes. As a result, we expect the frequency peak to be 
substantially broadened at such ambient temperatures. 

TABLE I: Calculated vibrational frequency shifts and inter- 
action energies for various orientations at the two types of 
binding site. The largest binding Channel A sites are at the 
intersections of that channel with the cell body diagonals (see 
Fig. [2]), while those for Channel B are at the cell boundaries 
(see Fig. [3]). The frequency shifts are relative to our calculated 
values for free H2 of 4160 cm~^. 
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Infrared absorption spectroscopy is particularly useful 
for studying the incorporation of H2 molecules in semi- 
conductors, as was demonstrated for amorphous silicon 
. The onset of IR activity and the position of the II2 
internal vibrational mode are both sensitive measures of 
H2 interaction with the matrix. The actual measure- 
ments are difficult because the absorption cross section 
of II2 is expected to be so weak that overtones and com- 
bination bands of the MOF material will interfere with 
the measurement. 

For these experiments, Zn2(BDC)2(TED) powder sam- 
ples were prepared solvothermally [j]. A small amount 
(10 mg) was lightly pressed on a KBr pellet, the consti- 
tutive solvent removed by activation (annealing for 10 h 
under 10~^ Torr vacuum) and the spectrum was recorded 
(using a clean KBr pellet as reference). The top panel 
of Fig. |4] shows the spectrum associated with overtones 
and combination bands of the clean MOF. Upon loading 
with high pressure He gas, these bands are perturbed, as 
evidenced by spectral shifts, due to minor rearrangement 
of the organic ligands. Consequently, in order to obtain a 
clean spectrum of II2 molecular absorption, D2 molecules 
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were used as a background for the H2 absorption exper- 
iments, as shown in the middle spectra in Fig. |4l The 
difference obtained is shown in the bottom three spectra 
for selected pressures (600, 700 and 800 psi). 

The main observation is the presence of a well-defined 
band centered at 4120 cm^^, which is composed of 3/4 
ortho and 1/4 para hydrogen (room temperature distri- 
bution), hence the asymmetry towards higher frequen- 
cies. The measured shift of —35 cm~^ from the unper- 
turbed ortho (4155 cm~^) and para (4161 cm~^) frequen- 
cies is in good agreement with the calculated —28/25 
cm~^ shifts. Furthermore, the integrated areas (mea- 
sured from 300 to 800 psi) are linear with pressure, which 
is expected since, even at 800 psi, the average H2 loading 
per unit cell is only '^l H2 out of ^^12 possible sites, well 
within the linear regime for sequential loading found by 
the calculations. The observed vibrational band is ~54 
cm^^ broad, consistent with phonon and thermal broad- 
ening at room temperature, and with tails extending over 
50 cm^^, which may arise from the translational states. 
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FIG. 4: Infrared absorbance of the Zn2(BDC)2(TED) MOF 
at room temperature. Top panel: clean MOF at ambient 
pressure; Middle panel: MOF at 800 psi of H2 (upper curve) 
or D2 (lower curve); Bottom panel: difference between the H2 
and D2 spectra at several pressures. 

To summarize, we have shown that a combination of 
experimental and theoretical methods gives a consistent 
and accurate picture of H2 binding in a prototype MOF 
structure. This picture will become still more precise 
when both low temperature infrared absorption can be 
performed, and the calculations of quantum effects on the 
motion of H2 within the MOF are complete. The meth- 
ods are now ripe for application to larger MOF structures 
with larger H2 binding energies and hopefully will give 



additional clues about structural and chemical changes 
that can increase these energies even further. 
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